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We show that the spin-orbit torque induced magnetization switching in nanomagnets presenting 
Dzyaloshinskii-Moriya (DMI) interaction is governed by a chiral domain nucleation at the edges. The 
nucleation is induced by the DMI and the applied in-plane magnetic field followed by domain wall 
propagation. Our micromagnetic simulations show that the DC switching current can be defined 
as the edge nucleation current, which decreases strongly with increasing amplitude of the DMI. 

This description allows us to build a simple analytical model to quantitatively predict the switching 
current. We find that domain nucleation occurs down to a lateral size of 15 nm, defined by the 
length-scale of the DMI, beyond which the reversal mechanism approaches a macrospin behavior. 

The switching is deterministic and bipolar 

PACS numbers: 75.60.Jk,85.70.Ay,85.75.Dd 


The recent discovery that a current can switch the 
magnetization of a nanomagnet in ultrathin heavy metal 
(HM)/Ferromagnetic (FM) multilayers has opened a new 
path to manipulate magnetization at the nanoscale [ij. 
The switching arises from structural inversion asymme¬ 
try and high spin-coupling, resulting in a spin current 
from the HM into the FM. This novel switching mech¬ 
anism has led to innovative magnetic memory concept, 
namely the spin-orbit torque MRAM [H-Q, which com¬ 
bines large endurance, low power, and fast switching and 
thus appears as a possible non-volatile alternative for 
cache memory applications. Recently, Garello et al. 0 
demonstrated deterministic magnetization switching by 
spin-orbit torque (SOT) in ultrathin Pt/Co/AlOa;, as 
fast as 180 ps. These observations could not be ex¬ 
plained within a simple macrospin approach, suggesting 
a magnetization reversal mechanism by domain nucle¬ 
ation and domain wall (DW) propagation. The failure 
of the macrospin approach for quantitative description is 
also underlined by the predicted switching current den¬ 
sity, which is nearly one order of magnitude larger than 
experimental ones [5l-l7|. Besides its fundamental impor¬ 
tance, this lack of a proper quantitative modeling is an 
important issue for the design of logic and memory de¬ 
vices based on SOT switchin g, which have so far consid¬ 
ered a macrospin description [8l-lll| . The final ingredient 
is the presence of antisymmetric exchange interaction, i.e. 
Dzyaloshinskii-Moriya interaction (DMI). This exchange 
tends to form states of non-collinear magnetization, pro¬ 
moting homochiral Neel DW fl2l - [l3| . In the Neel config- 
uration, a maximal SOT is applied on the DW [1^ IlSl - 


17| . This explains the large current induced DW velocity 


observed experimentally tllllTl. Moreover, the DMI can 
result in significant magnetization tilting at the edges of 
magnetic structures, resulting, e.g., in asymmetric field 


induced domain nucleation 18. JH. The influence of the 
DMI on the magnetization pattern during SOT switch¬ 
ing was recently pointed out in micromagnetic simulation 
studies [20l - [^ . whereas recent experimental work 231 
explained the SOT switching mechanism by the expan¬ 
sion of a magnetic bubble. 


Here, using micromagnetic simulations and analytic 
modeling, we show that the SOT-induced magnetization 
switching in the presence of DMI is governed by domain 
nucleation on one edge followed by propagation to the 
opposite edge. This reversal process allows to explain 
the ultra-fast deterministic switching observed experi¬ 
mentally. We systematically demonstrate that DMI leads 
to a large decrease of the switching current and of the 
switching time and thus strongly affects the reversal en¬ 
ergy. On the basis of our micromagnetic simulations, we 
provide a simple analytical model which allows to quan¬ 
titatively predict the SOT switching current in the pres¬ 
ence of DMI. Finally, we address the evolution of the 
switching mechanism as the lateral dimension decreases, 
which is a key feature for the device scalability. 


The structures considered in this study are similar to 
the one used by Garello et al. Q : a perpendicularly mag¬ 
netized Co circular nanodot on top of a Pt stripe and 
capped with alumina. The DMI is included into the simu¬ 
lation using the expression of Ref. [l^ . In addition to the 
standard micromagnetic energy density (which includes 
the exchange, the magnetocrystalline anisotropy, the Zee- 
man and the demagnetizing energy), the current injected 
in the Pt layer leads to two SOT terms in the the Landau- 
Lifshitz-Gilbert equation: the field-like Tfl oc m x and 
the damping-like Tdl oc to x (to x e^), where Cy is the 
unit vector in y-direction (see 2^ for additional details). 
If not mentioned otherwise the external applied field is 
PoHiipp = —0.1 T, and the material parameters are (^ : 
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FIG. 1. (color online) a) Time evolution of the average out-of-plane magnetization for different applied current densities 
(variations in steps of 10^^ A/m^). The minimum current to trigger switching, i.e. the critical current, is highlighted in blue. 
The green curve indicates the threshold of turbulent behavior (see text), b) Sketch of the magnetization conhguration at 
different stages of the switching process, c) Magnetization orientation in the center and at the left and right edges. The current 
induced damping-like torque (represented as effective field IFdl) can only drive the left edge magnetization into instability, 
resulting in a nucleation at the left edge, d) Snap-shots of the magnetization conhguration showing the reversal from down 
(black) to up (white) via domain wall nucleation and propagation under an externally applied held of fioH = 0.1 T and a 
current density of 2.6 10^^ A/m^. 


the saturation magnetization Ms = 1090 kA/m, the ex¬ 
change constant Aex = 1-0 x 10“^^ A/m, the perpendicu¬ 
lar magnetic anisotropy constant = 1248 kJ/m^, the 
DMI amplitude D = 2 mJ/m^, the Gilbert damping pa¬ 
rameter a = 0.5, the torques TSj = —0.05 pTm^/A and 
= +0.1 pTmVA. 

The 3D micromagnetic simulations are performed us¬ 
ing the solver MicroSD with a mesh size smaller than 
1.5 nm. The initial magnetization state of the dot is the 
remanent state after saturation by a negative magnetic 
field (—Oz) as shown for 0 ps in Fig. [T]d). In the pres¬ 
ence of an applied magnetic field i?app in the x-direction, 
magnetization dynamics is induced by a current pulse 
with a rise (and fall) time of 50 ps and variable width 
and amplitude. Typical simulation results of a 100 nm 
dot are presented in Fig.[l]a). Depending on the current 
amplitude, three regimes are identified: 

1) For Japp < 2.50 10^^ A/m^ no magnetization switch¬ 
ing is observed. The SOT leads to a slight tilting of 
the magnetization toward the plane of the dot, but the 
magnetization relaxes toward its initial equilibrium state 
after the pulse. 

2) At intermediate current values (2.60 10^^ A/m^ < 
Japp < 3.70 10^^ A/m^) magnetization reversal occurs. 
The time evolution of the magnetization pattern in the 
dot, see Fi g.fj d) reveals that, in contrast to recent inter¬ 
pretations [231 ) , the magnetization reversal occurs by do¬ 
main nucleation shortly after the pulse injection (100 ps), 
followed by fast DW propagation. The nucleation always 
occurs on the left edge of the dot. Once nucleated, the 


DW propagates fast through the dot and is expelled on 
the opposed edge. The switching time to, defined by 
(mz){t = to) = 0, decreases as Japp increases; the in¬ 
crease of the slope of {mz){t) indicates that this is re¬ 
lated to a faster DW propagation. As expected the DW 
has a Neel configuration due to the large DMI. The sim¬ 
ulation highlights that the DW nucleation occurs for all 
current values on the same edge in a deterministic way. 
Symmetrically, when reversing the sign of the current, 
the reversal from the up to the down state occurs on the 
opposite edge, i.e. the behavior is bipolar. 

3) For higher currents (Japp > 3.70 10^^ A/m^) the 
motion of the DW becomes turbulent (oscillatory) and 
the coherence of the switching is destroyed. 

The magnetization reversal scheme can be explained in 
a simplified manner by considering the combined effect of 
DMI, external magnetic field, and SOT, but neglecting 
small variations of the demagnetizing field [i^ . The DMI 
is too small to introduce a spin spiral but results in a 
magnetization canting at the dot edges [Hill,nil. The 
edge canting can be considered as an additional effective 
field with spacial variation: on one side this field adds 
to the in-plane applied field, while it counteracts on the 
other (see Fig. fflb). This leads to an asymmetric tilting 
of the magnetization on both edges. 

Upon current injection the damping-like torque 
emerges. Its effect can be interpreted as a rotating mag¬ 
netic field of the form T/dl oc Jappiri x e/ (see Fig. [T]c)). 
This leads to a rotation of the magnetization towards the 
film plane on one side and away from the film plane on 
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FIG. 2. (color online) a) Critical current for destabilizing the system as a function of DMI strength, b) The relation between 
the critical current and the switching time to for two different values of DMI.The inset show the data for D — 2 mj/m^, but in 
a transformed coordinate system Japp verus . c) The switching time versus current for different dot diameters, d) Critical 
current and current density for different dot sizes. The calculation of the current assumes a 3 nm thick Pt line. 


the other. Naturally, the current polarity is chosen such 
that the stronger tilted edge magnetization turns towards 
the film plane. Above a critical current an instability oc¬ 
curs, leading to domain nucleation and consecutive DW 
propagation. It is clear that the current required to 
introduce the instability, reduces with increasing DMI. 
This behavior is seen in Figl2]a), where Jc tends to zero 
when iJ ss 3.8 mJ/m^. Moreover, for Japp > Jc an in¬ 
crease of DMI decreases the switching time, as can be 
seen from Fig. [2]b). After expelling the DW on the op¬ 
posite side, switching has occurred and the more tilted 
edge appears on the opposite side. As the SOT rotates 
this side away from the film plane and is not sufficient 
to rotate the less tilted side into instability, the state is, 
hence, stable. It can be easily checked that this rever¬ 
sal scheme is in agreement with the hysteretic bipolar 
switching observed experimentally when sweeping Japp 
and iJapp 0 • 

To understand these results better, we consider a sim¬ 
ple analytical model, which describes the reversal pro¬ 
cess in the presence of both DMI and SOT. Using a La- 
grangian approach and following Pizzini et al. [l9|, the 
strategy is, eventually, similar to the Stoner-Wohlfarth 
approach in a single domain particle but using the en¬ 
ergy functional 


m 

V 


= — ATeffCOS^ 0-MsiJappSin0 —A/sJappIbL^*, (1) 


where the effect of the SOT is introduced by the last 
The equilibrium magnetization angle in 


24|. 


term 

the center, 9c is found by minimizing Eq. 1161 while 
edge angle 9e is found by solving [E{6e) — E{9c)]/V = 
0.25D^A~^ [ 2 ^. For small SOT and Happ, two stable so¬ 
lutions for 9e exist, corresponding to both sample edges. 
Above a threshold SOT one solution disappears, indi¬ 
cating that the magnetization on one edge is unstable, 
i.e. domain nucleation occurs. Using numerical meth¬ 
ods, the critical current for nucleation Jc can be cal¬ 
culated easily as a function of D (see Fig. [5] a), black 


line). A good agreement is obtained with micromag- 
netic simulation for a dot diameter d = 100 nm (circles). 
For D tending to zero, the nucleation current tends to 
the critical current predicted by the macrospin model 
Jc = 4.1 X 10^^ A/m^ Q. The absence of full quanti¬ 
tative agreement with micromagnetic simulation can be 
attributed to variations of the demagnetizing tensor and 
variations of the magnetization along the y-direction due 
to the curvature of the dot. Better agreement is obtained 
when neglecting these effects in a quasi ID simulation 
(square dots). Note that this nucleation current is actu¬ 
ally the threshold current for quasi-DC current pulse. 


In the following, we discuss the dynamics of the mag¬ 
netization switching. In Fig. [JJd) the switching time is 
shown as a function of Japp > Jc- With increasing Japp 
the switching time decreases rapidly as the DW velocity 
increases [1^. If D is reduced, the DW propagation is 


slower, resulting in a larger switching time. In the inset 
we show Japp versus 1/to for D = 2 mJ/m : a linear 
scaling is observed, in qualitative agreement with exper¬ 
iment yi. 


Naturally, to depends on the dot diameter. This is 
a key parameter for SOT applications. The evolution of 
switching time versus current density for varying dot sizes 
is shown in Fig.[^c). When decreasing the diameter from 
100 nm down to 50 nm, a shift to shorter switching times 
is observed while a slightly higher onset current is found. 
Similar behavior is found when decreasing the size down 
to 30 nm and further down to 25 nm. It is, however, im¬ 
portant to note that the latter two graphs become iden¬ 
tical for larger Japp. Reducing the size down to 15 nm, 
results in a dramatic increase of the threshold current 
density. Moreover, deterministic switching is observed in 
a narrow current density region only. Overall one has 
indications for three different size-dependent switching 
regimes. In the first regime the switching is covered by 
nucleation and propagation of a DW and the decrease of 
to is mainly caused by a reduced distance for the DW 
to travel. In the second regime the switching remains 
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governed by DW propagation. The diameter, however, 
becomes comparable to approximately twice the value of 
^ = 2AID K, 10 nm, the characteristic length scale on 
which canting of the edge magnetization is observed. In 
this situation the edge angle due to DMI differs from 
the ideal infinite case and opposite ed ges are not com¬ 
pletely independent anymore (see Ref. 2^). While this 
does not cause coherent rotation yet, it affects the DW 
motion. The coherent regime is reached at diameters in 
the range of the DW width A = k, 14 nm. 

This explains the significant change in switching behav¬ 
ior for the 15 nm dot. Note that the switching current at 
this size is close to the one predicted by macrospin sim¬ 
ulation (4.1 X 10^^ A/m^). It is worth mentioning that 
while the current density strongly increases with decreas¬ 
ing dot diameter, the current in the 3 nm thick Pt stripe 
decreases almost linearly, as can be seen from Fig. [2]d). 
Therefore, the device exhibits favorable scaling behavior 
and assuming a 1 kD resistance for the addressing tran¬ 
sistor of a 30 nm dot, switching in about 300 ps, needs 
only 20 fJ for one switching event, which is significantly 
smaller than the energy for perpendicular spin-transfer 
torque devices [28 1. 

Naturally, the threshold current and switching time 
depend on several intrinsic as well as extrinsic parame¬ 
ters. We have studied in detail the influence of the ap¬ 
plied field, the damping constant, the strength of the field 
like torque, and temperature. The results are shown in 
Fig. [31 Variations in these parameters lead to quanti¬ 
tative changes of the nucleation current as well as the 
switching time. In all cases this is mainly attributed 
to changes in DW velocity; lower damping increases the 
wall velocity and so does an in-plane field, as it promotes 
and stabilizes a Neel type wall. A negative field like 
torque also stabilizes the DW, while a positive one desta¬ 
bilizes it, therefore increasing the switching time. The 
edge nucleation/DW propagation mechanism, however, 
is not affected. Most importantly Fig.|3]a) shows that the 
mechanism of switching by nucleation and propagation is 
very robust against fluctuations due to temperature (See 
Ref. 2^ for more details). The temperature fluctuations 
strongly decrease the threshold current (Fig. |3|b). Tem¬ 
perature effectively lowers the nucleation barrier, such 
that nucleation times get shorter and, consequently, the 
whole switching becomes faster. It has to be pointed out 
that the nucleation still takes place at the same position 
on the dot edge and the overall process remains bipolar 
with respect to field and current reversal. This temper¬ 
ature robustness, however, strongly relies on the large 
damping, as can be seen from Fig.|3|b). With decreasing 
a an increasing tendency of oscillations is observed, such 
that deterministic switching cannot be guaranteed 

To conclude, we have studied the current induced 
magnetization switching of a nanomagnet by spin-orbit 
torques in the presence of Dzyaloshinskii-Moryia interac¬ 
tion (DMI). The critical switching current strongly de- 
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FIG. 3. (color online) a) Switching time as a fnnction of ap¬ 
plied current density, varying intrinsic and extrinsic param¬ 
eters. For the temperature case, the average to is plotted. 
The single event switching time is defined as before, while 
the average to is defined as the time where the probability 
of stochastical switching reaches 90%. b) Several switch¬ 
ing graphs {mz){t) for varying damping at T = 50 K and 
dapp = 2.6 10^^ A/m^. For fixed a variations are only due to 
temperature fluctuations. 


creases with increasing amplitude of DMI and we pro¬ 
vide a simple analytical model for this dependency. This 
switching mechanism via chiral domain nucleation ex¬ 
plains the deterministic switching observed experimen¬ 
tally in ultra thin Pt/Co/AlOa, even for sub-ns pulses. 
The switching is mainly introduced by the damping like 
torque, but the field like torque cannot be neglected as 
it strongly influences the switching time. Our systematic 
studies show a change in the reversal mechanism below 
diameters of 30 nm, while the switching remains deter¬ 
ministic and bipolar. However, at 0 K the operational 
window for current densities decreases with decreasing 
dot diameter. The influence of temperature on this tech¬ 
nological important limit will be investigated in the fu¬ 
ture. Most importantly, current scalability is maintained. 
Confirming the potential of SOT-MRAM for scalable fast 
non-volatile memory application, our results will help in 
the design of devices based on this technology. 
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Appendix 

General description of sample and technics 


energy of the form 

Me; • 

Ekin =- -(pcosO, (4) 

7 

where the dot-notation refers to the total time derivative. 
With U as potential energy, the equation of motion is 
derived from the action 


The properties of cobalt films sandwiched between 
platinum and a insulating oxide such as AlO are stud¬ 
ied. Due to the different substrate and capping mate¬ 
rials the inversion symmetry is broken along the ver¬ 
tical axis (Oz). The magnetization is oriented out-of- 
plane with a strong magnetocrystalline perpendicular 
anisotropy. In addition to the standard micromagnetic 
energy density, which includes the exchange, the magne¬ 
tocrystalline anisotropy, the Zeeman and the demagne¬ 
tizing energy, the Dzyaloshinskii-Moriya contribution is 
included according to the following relation: 


Edm = D 


m. 


dmx 

dx 


drriz 

dx 


dniy 

dy 


dm. 


dy 


( 2 ) 

This expression corresponds to a sample isotropic in the 
x-y-plane, where the Dzyaloshinskii vector for any in¬ 
plane direction is D{ez x Cu) with D a uniform con¬ 
stant, originating from the symmetry breaking at the z- 
surface 13|, and being the unit vector in direction of 
an arbitrary u. 

Micromagnetic simulations are based on the time inte¬ 
gration of Landau-Lifshitz-Gilbert equation including the 
field-like Tfl and damping-like Tdl spin-orbit torques: 


dm 

dt 


= -7o 


m X + Hdm) 


-|- a m X 


dm 

dt 


+7fl + Tdl 


(3) 


Here m is the unitary vector of the magnetization, a 
is the Gilbert damping parameter, 70 is the product 
of the vacuum permeability yo and the free electron 
gyromagnetic ratio 7 . For the spin-orbit torques only 
the first order terms are considered, namely Tfl = 
7T^L'^app(hT. X Cy) and Tdl = yT^L'^apphi x (m x e^), 
where Tp^ = —0.05 pTm^/A and T^p = -1-0.1 pTm^/A 
are scalar constants. Higher order terms can be found in 
Garello et al. [j^ . 

The appearance of Cy in the definition of the torques 
is a simplification due to currents in Cx- In the general 
case this must be replaced by Cy ^ Cu = Sz x e), where 
Bj is the unit vector in the direction of the conventional 
current, i.e. opposite to the electron flow. 


Lagrange formalism for damping like spin-torque 


E = Tkin — G, 


(5) 


in the typical form of 

^ _ _d ^ ^ 

6q df Sq Sq ’ 

where 6 refers to the functional derivative and q € {0, </)}. 
Here we add a dissipative term F of the form 


aMs 


— m -I —m X 
\dt a 



(7) 


and X = 77DL'7app- The square of the hrst term in paren¬ 
thesis of Eq. [3 results in the standard damping of the 
LEG equation. The square of the second term vanishes 
when evaluating Eq. [ 6 l while the mixed term results in 
the damping-like torque. Using a quasi ID case with 
potential energy of the form 

U = A6l- DOx - Aeffcos^ e - yoMsH^pp sinB. (8) 


(here the index x indicates the partial derivative with 
respect to x) and a field applied in a:-direction one derives 
the equations of motion as. 


Ms • 

- ipsinO — Kessm29 -|- 2A6xx+ 

7 

+yoMsHs,ppCos9 -^0-1-^xcos0 = O (9) 

7 7 


Ms 


9 sin 9 


aMs 


(j) sim 0 - 1 - 


Ms 


X sin (j) sin 9 cos 0 = 0 (10) 


Note that the DMI drops out when applying the func¬ 
tional derivative. In the quasi static case oi 9 = (j) = 0 
this simplifies to 


- A:effsin20 -I- 2A9xx + yoMsH^pp cos 9+ 

+ XCOS()> = 0 ( 11 ) 

7 

Ms 

-X sin (j) sin 0 cos 0 = 0 (12) 

7 

The second equations is already consistently fulfilled by 
setting (j) = mr, n G Zq, i.e. the quasi ID case in the 
x-z-plane. It remains to solve 


The LEG equation in spherical coordinates (0, (j)) can 
be derived from a Eagrangian by defining a pseudo kinetic 


— KeS sin 20 -|- 2A9xx + /roM^giEapp 


Ms 

cos 0 ±-x = 0 (13) 
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FIG. 4. a) Angle of the magnetization with respect to the «-direction as a function of an applied magnetic field in ^-direction 
and and zero current. Due to DMI the magnetization at the edges of the infinite long stripe differ from the center. At 
approximately 0.5 T the first edge is fully in-plane. Note, without current and due to symmetry the maximum edge angle 
is always ti/2. b) Analytically calculated critical current Jc, required to destabilize the system, plotted against an in-plane 
applied magnetic field. In both cases the analytical model (continuous lines) is well reproduced by a micromagnetic simulation 
(dots) for a stripe of width 200 nm (128 cells) and a constant effective anisotropy. The other parameters are given in the text. 


Multiplying by Ox allows to integrate resulting in 

AOl = -Keff cos^ 0 - /roMsiJapp sin 0 =F —+ C, (14) 

7 

where C is the integration constant. This result is 
identical to Ref. [l^ except for the linear term due to 
the damping like torque, i.e. g{Japp)0 with g{Japp) = 
-^S'^app^DL- 


Analytical model for the critical current at 0 K 


A simplified analytical model has been developed, 
which already exhibits the important mechanisms of the 
behavior observed in experiment. An important ingredi¬ 
ent for deterministic switching in the analytical model is 
the presence of DMI. The DMI in the presented material 
system is below the critical value, such that no spin spiral 
is formed. Due to the boundary condition 


D ^ ^ 

— {Cz X n) X m = Vftm 


(15) 


and the finite size of the sample, however, the magneti¬ 
zation is not homogeneous and canting is observed at the 
edges. Here n is the surface normal and Vfj the normal 
derivative. It has been shown recently that the canting 
angle at the edge as well as in the center of the sam¬ 
ple can be evaluated using simple energy arguments 0, 
provided that the demagnetizing energy can be approx¬ 
imated by an effective anisotropy. This is not always 
the case [27|, but for the given parameters it is justified 
as can be seen from Rohart and Thiaville 18|. The de¬ 
rived angles from this simple analytical calculation and a 
micromagnetic simulation that accounts for the demag¬ 
netization energy by an effective anisotropy, are in very 


good agreement, as can be seen from Fig. [4] a). As the 
micromagnetic angle is taken from the cell magnetiza¬ 
tion, a small deviation to the theoretical edge value is 
observed. With decreasing mesh size this deviation de¬ 
creases. The original publication 0 calculates a critical 
magnetic field in ^-direction resulting in the nucleation of 
a domain. In the present case the according term is sub¬ 
stituted by a current dependent term (see paragraph p), 
taking account for the dominant damping like torque. At 
first glance this is problematic, as the model considers a 
static case approaching an instability. Moreover, it is 2D 
in the sense that it assumes my = 0, while a damping 
like torque initially acts along the y-direction. For large 
damping, however, changes of the magnetization become 
quasi static and mainly take place in the a;-z-plane. 

To provide the full solution for 0 = 0(x) one needs 
to integrate Eq. (TD In the center of the sample, i.e. 
far away from the edges, one can assume Ox = 0, such 
that the DMI drops out again. The total energy density, 
hence, has the form 

= -RieffCos^ 0 - MsHapp sin6» - g{Japp)0 (16) 

and the equilibrium angle is found by minimizing Eq. 1161 
A second position where one can find a solution is the 
edge. Here additional information is given by the bound¬ 
ary condition Eq. 1151 Inserting the boundary conditions 
into Eg. [HI results in Eg. [HI plus an additional offset AE 
of 


AE _ 

~ ~ 4A 


(17) 


As a result the solution of the edge angle is given by 
a value 0 that has a AE higher energy than the mini¬ 
mum energy in the DMI free energy landscape of Eq. 1161 
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FIG. 5. Accumulated switching probability as function of 
time for different current densities at T = 50 K. 

Hence, a stable solution for the edge angle can only be 
found if the energy landscape provides a value of /S.E 
above the local minimum, i.e. above the solution of the 
center angle. For a small current the edge solution cor¬ 
responds to a point near the next local maximum. With 
increasing current this local maximum decreases and at 
the critical current Jc, it is separated from the minimum 
by exactly AE’. Eventually, the solution for the mag¬ 
netization angle at the sample edge vanishes for larger 
currents, i.e. the current drives the system into instabil¬ 
ity. The results of the analytical model for current driven 
instabilities are plotted as continuous line in Figure [4] b). 
The analytical values are compared to a simple simula¬ 
tion neglecting variations of the demagnetizing field and 
assuming only an effective anisotropy. The agreement of 
analytical results and simulation are in astonishing agree¬ 
ment. The applicability of the model is, however, limited 
as it neglects temperature. It is shown below that tem¬ 
perature has a significant effect on the critical current. 



FIG. 6. Cross section showing for different dot diameters 
and a 100 mT applied field x-direction. Within approximately 
20 nm the canting at the edge decays towards a homogeneous 
magnetization at the dot center. 

the integrated probability reaches 90%. Naturally, the 
switching time decreases with increasing current density. 
One has to keep in mind, however, that this data rep¬ 
resentation does not give information about oscillatory 
behavior for larger times beyond the (first) zero crossing 
of {ruz)- This graph, hence, can pretend fast switching 
where actually oscillatory behavior is present. 


Thermal fluctuations and switching probabilities 

The temperature was included in the form of a Gaus¬ 
sian distributed thermal field 77th: which is added to the 
effective field HeS- The thermal fluctuations have the 
following properties [s^: 

= 0 ( 18 ) 

and 

(i7th.*(t)i7th.,(t')) = (7oMs(t - t')) (19) 

Mo Fcell 

where /cb is the Boltzmann constant and Well the vol¬ 
ume of the discretization cell. As temperature results 
in stochastic behavior of switching one has to consider 
switching probabilities. This is done by evaluating 100 
independent switching events for each parameter set. For 
each event the the switching time to, i.e. the time when 
{niz) crosses zero is determined. From this the integrated 
probability is calculated. A typical result is shown in 
Fig. m The overall switching time in the presence of 
thermal fluctuations is then defined as the time where 


Size dependence 

In addition to the exchange length, the system at hand 
exhibits a second important length scale. This length is 
given by the interplay of DMI and exchange interaction, 
i.e., to what extend the non-collinear magnetism from 
the edge penetrates the sample. The according length is 
given by ^ = 2A/D « 10 nm. It can be seen from Fig. [6] 
that this length scale becomes relevant for dot diameters 
below approximately 40 nm. For smaller diameters a 
strictly monotonic change of rUz is observed while larger 
diameters present a center plateau. 
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